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Abstract 

The present article reveals important properties of the confluent Heun's functions. We derive a set of 
novel relations for confluent Heun's functions and their derivatives of arbitrary order. Speciflc new subclasses 
of confluent Heun's functions are introduced and studied. A new alternative derivation of confluent Heun's 
polynomials is presented. 

PACS: 02.30.Gp, 02.30.Hq 



1 Introduction 

The solutions of the confluent Heun's differential equation fTI-IB], written in the simplest uniform shape [7] 

H''+(a + ^ + l±l)H'+(^ + ^)H = 0, (1.1) 
\ z z — 1/ \z z — 1/ 

are of continuous and significant interest for many appUcations in different areas of natural sciences and especially 
in physics. 

The equation (|l.ip has three singular points: two regular ones - z = and 2 = 1, and one irregular one - 
z = oo. The standard confluent Heun's function HeunC(a, /3, 7, i5, 77, z) is a unique particular solution, which is 
regular around the regular singular point z = 0. It is defined via the convergent in the disk |z| < 1 Taylor series 
expansion 

00 

HeunC(a,/3,7, (5, ?7, z) = ^ w„(a, /3, 7, (5, 77)2", (1.2) 

ri=0 

assuming the normalization HeunC(a, /3, 7, (5, 77, 0) = 1. The parameters a, /?, 7, (5, ry, introduced in [31 13] and 
used in the widespread computer package Maple, are related with fi and v according to the equations /i = 
^{a - j3 - J + a(3 - Pj) - 1] and ly ^ ^{a + (3 + j + + Pj) + S + ?]. 

The coefficients w„(a, /3, 7, (5, 77) are determined by three-term recurrence relation 

AnVn = BnVn-1 + C„W„_2 (1.3) 

with initial condition v_i =0, vq — 1. Here 

An = 1 + — — > 1, when n oo, 

77- 

, , -a + /3 + 7-l , 77 -(-Q + /3 + 7)/2- a/3/2 + /37/2 

i>„ = 1 H 1 1, when ^ oo, 

77 

a / S /9 + T \ 

C„ = [ - + — TT^ + n - 1 ^0, when ti ^ oo. (1.4) 
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According to [T]-[6] the function HcunC(a, /3, 7, S, rj, z) reduces to a polynomial of degree > with respect 
to the variable z if and only if the following two conditions are satisfied: 

S /3 + 7 

- + -—^ + N + 1^0, 1.5a 
a 2 

AN+i{fi)=0. (1.5b) 

Further on we call the first condition (|1.5a[) - a " JAr-condition" , and the second one (jl.5bp - a " AAr+i-condition" . 
One can find an explicit form of the left hand side Ajv+i(m) of the condition (|1.5bp . convenient for practical 
calculations, in the Appendix. 

Indeed, the (Sjv-condition is equivalent to the equation Cn^2 = 0; f^nd the Ajv+i -condition turns to be 
equivalent to the requirement ujv+i (a, (3,j,S,r]) = 0. Then as a result of equation p.3p and additional conditions 
(jl.Sp we obtain ujv+2(a7 /?, 7, (5, ry) = 0. Since two consecutive terms in the tree-term recurrence relation (|1.3p 
are zero, all next terms are zero, too. Hence, under simultaneous fulfillment of the two additional conditions 
(jl.Sp the confluent Heun function HeunC(Q;, f3, 7, S, rj, z) (|1.2|1 reduces to a polynomial of degree A^. 

In present article we derive a set of novel relations and differential equations for confluent Heun's functions 
p.2p and their derivatives ^^HeunC(a, 7, (5, 77, z) of arbitrary order n - see Section 2. In Section 3 we 
introduce a new subclass of confluent Heun's functions HeunCAr(Q;, /?, 7, 77, z), which obey only the (^Ar-condition 
(|1.5ap . as well as their associate confluent Heun's functions HeunC'^(a, /?, 7, 77, z). In Section 4 we utilize the 
newly found relations to present an alternative derivation of the confluent Heun's polynomials without use of 
the recurrence relation (|1.3p . Thus our consideration reveals important new properties of the confluent Heun 
functions. 

In the concluding Section 5 we briefly discuss the relation of our results for confluent Heun's functions with 
similar results for general Heun's functions. 

The confluent Heun's functions HeunC(a, /3, 7, 5, 77, z), the specific classes of i5jv-confluent Heun's functions 
HeunCAr(Q;, 7, 77, z) and their associated functions HeunC^(Q;, 7, 77, z), as well as confluent Heun's polyno- 
mials PHeunCAr_fc(a, /3, 7, z) play a very important role in some applications, especially in gravitational physics 
[7]-[16j. The present study of their properties was inspired by our desire to reach a true mathematical and 
physical understanding of the important Teukolsky-Starobinsky identities, derived and used in the early arti- 
cles [1Z]-[19]. We present here the formal mathematical results, which can be applied, too, in other scientific 
domains, both for analytical and for numerical calculations. 

2 Novel relations for the confluent Heun's functions and their 
derivatives 

Let us define the differential expression 

D^,0,,,s,, ^z{z-l)^—+{^a+^ + —j) _ + - + __ j . 
One can use it to write down the confluent Heun equation (jl.ip in the following compact form 

It can be shown that certain restrictions of differential expression -Da, /3. 7, 5, 77 on proper functional spaces yield 
self-adjoint differential operators [5]. In the present article we will skip the detail, needed for justiflcation of the 
operators' domains and proper scalar products in the corresponding linear spaces of functions. Here we restrict 
our consideration only to formal manipulations with differential expressions like (j2.ip . 

The confluent Heun's operator Da,f3,j,5,ri (|2.ip owns a remarkable property. Its eigenfunctions Hx{z) for 
eigenvalues A 7^ 0, i.e. the solutions of the ordinary differential equation 

Dc.,i3,jAriHx{z) = XH),{z) (2.3) 
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(2.1) 
(2.2) 



are at the same time solutions of confluent Hcun equation (|2.2[) with the same parameters a, (3, 7, S and a 
different parameter rj' = 7] — \. For them the following confluent Heun's equation takes place 

xHx{z) = 0, VA e C. (2.4) 



Indeed, the equation (j2.3|) obviously can be rewritten in the form (|2.2p with fi' ~ fi + X and v' = i> — X 
ten, 
result. 



Then, using the relations S = i> + fi — a — 1^ and rj — — /3 — 7 + a(3 — (3"f) — fi one obtains the above 



Commuting Da,fj.-y,5,ri with differential expression we derive the basic novel relation 

d" ~ - d" /J /? + 7 \ <i"-i 

-^Da.,(},^,5,r, = £'a(„),/3(«),7(«),<5(n),,,(«)^ + ntt I - + — ^ + n 1 ^^^;irT ' n = 0, 1, 2, . . . (2.5) 
The transformation {a, /?, 7, (5, 77} — + {a(n), /3(n), 7(n), (5(n), ?7(?t.)} defined by the relations 

(^(71) (5 Ti 

a{n) = a, /3(n) = /3 + n, 7(n) = 7 + n, — — — - = V n, r]{n) = rj + —{n — a + (3 + j) (2-6) 

a{n) a 2 

yields a specific augmentation of the indexes a, (3, 7, S and rj. 

To simplify our notations we will use further on the following ones: 

£'q,/3,7,<5,j) = Do, £'q(1),/3(1),7(1),5(1),j7(1) = • • ■ , Da(n),f3{n),j{n),5(n).ri{n) = Dn, ■ ■ ■ (2.7) 

Applying the equation (|2.5p to arbitrary solution i?(z) of confluent Heun's equation (|2.2p we obtain the 
following four-terms recurrence relation for the derivatives ■j-^H{z) — iJ'^"^(z): 



L>„i/(")(z) = -na (^^ + ^^-^ + H^''-^\z), n = 0, 1,2, , 



(2.8) 



For 71 > 2 it can be considered, too, as an ordinary differential equation of third order for the functions 
i7("-i)(z). 

Then applying several times the equation (12. 8p we obtain the following sequence of relations for arbitrary 
solution H{z) of confluent Heun's equation (|2.2p : 

D^D^...br.H^^\z)^{-arn\(- + ^4;l + l\H{z), ri = l,2,... (2.9) 

Here (a:)„ = F ( a: + 7i)/r (a;) = x{x + 1) . . . (a; + 7i — 1) stands for the Pochhammer symbol [5]. 

Finally, applying Dq to the both sides of equation (|2.9p we end with the ordinary differential equations of 
order 2{n + 1) for derivatives H^^^\z) of the solution H{z) of confluent Heun's equation (|2.2p : 



DoDi---DnH'-''\z) = 0, 71 = 1,2,... (2.10) 

The relations (|2.9p obviously show that any solution H{z) of confluent Heun's equation (|2.2p is an eigen- 

function of the operator D1D2 ■ ■ ■ Dn^jr: and (—a)" ti! + + 1^ is the corresponding eigenvalue. The 

relations (|2.10p show that the solutions H{z) of the confluent Heun's equation belong simultaneously to the 
null-spaces of the infinite sequence of linear operators DqDi . . . Dnjpr'- n = 0, 1, 2 . . . 

3 New subclass of confluent Heun's functions 



When the (^Ar-condition (|1.5ap is fulfilled, one obtains S = —a yl^Y^ + N + lj = 6n for some fixed nonnegative 
integer > 0. Then for any function H{z) G the equation (|2.5p reduces to 



D,H\=DM+A-r-^\. (3.1) 
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This relation shows that in the case under consideration the operator ^^n+i defines a specific generalized 
Darboux transformation [20 for the confluent Heun's operator Dq. It is remarkable that in the right hand side 
of equation p.ip we have just another confluent Heun's operator - Z3jv+i- As a result, if the function H\{z) 
with (5 = (5 AT is an eigenfunction of the confluent Heun operator Dq for some eigenvalue A, then the relation 
(j3.ip shows that the function ^^n+i H\{z) is an eigenfunction to the corresponding operator Dn+i for the same 
eigenvalue A. 

Note that the condition 6 — 6n invariant under arbitrary changes of the eigenvalue A. The operator Dq 
with S = Sn preserves the linear envelope of its eigenf unctions with different A, with the same 6 = 5n and fixed 
nonnegative integer N. The discussion of the equation (|2.3p at the beginning of the previous Section 2 shows 
that studying this linear space of functions, further on we can restrict our consideration with the case A = 0. 
All formulae for A 7^ can be obtained from the ones for A = replacing rj with rj — X. 

Suppose the function H{z) is in addition a solution of the equation (|2.2p . Then the relation (|3.ip yields the 

following novel confluent Heun's equation for the derivative '^^^4.^ : 

/ jN+l TT\ 

i)..r(^)=0. (3.2) 

Consider the two unique solutions of the Heun's differential equations (|2.2p and (|3.2p , which are simultane- 
ously regular at the point z = and equated to unity at this point. Then these solutions are precisely the two 
confluent Heun's functions HcunC(a, /?, 7, 6n, 77, z) and HcunC(a(A^+l), l3{N+l), j{N+l), Sn{N+1), r]{N+l), z), 
the last being constructed according to equations ()2.6[K For them the uniqueness theorem and the uniform nor- 
malization at the point z = yield the novel relation |^ 

dzN+i HeunC(o;, /3, 7, Sn, ??, z) = (TV + 1)! Ujv+i(a, 7, (^w, ??) x 

X HeunC(a(iV+ 1), f3{N+ 1), j{N+ 1), Sn{N+ 1), t]{N+ 1), z). (3.3) 



Here we have used the (N + l)-th coefficient in the Taylor series expansion (|1.2p . 

Definition 1: We call 6];^ -confluent- Heun- functions the functions HeunC(a, /3, 7, (Sat, 77, 2:) which obey the 
^AT-condition for some specific nonnegative integer N > and denote them as 

HeunCAr(Q;,/3,7,77, z) = HeunC(a, ^9, 7, -a{N (/3 + 7)/2) , 77, z). (3.4) 

Definition 2: Wc call associate (^Ar-confluent-Heun's-functions the functions 

HeunC* (a, /3, 7, 77, z) = (3.5) 

^}ieunC{a,(3+N-hl,-/+N + l,-a{(3+j)/2,f]-h{N+l)(N+l-a+(3+i)/2,z). (3.6) 

Now the relation p.3p can be represented in the short form 

^^jv+i HeunCAr(a, /3, 7, 7?, z) = Ta^ (a, /3, 7, 7?) HeunC* (a, /3, 7, 7?, z), (3.7) 

where we are using the specific constant 

J'jv(a,/3,7,^) - (iV + l)!«Ar+i(a,/3,7,-a(7V+l + (/3 + 7)/2),77). (3.8) 

According to the relations ((2^ we obtain S^iN + l)/a{N + 1) + {/3{N + 1) + 7(A^ + l))/2 iV + 1 > 0. 
Thus the (5jv-condition is not fulfilled for the associated J^v-confluent-Heun's- function HeunC* (a, /?, 7, r;, z). 
Hence, it does not belong to the class of the (5jv-confluent-Heun-functions. 

^The explicit form the relation 1)3. 3|l reads 

^^Kcunc(^a,l3,-y,-a(^^ + N + l^ ,7j,z^ = (N + iy.VN+i(a, (^^ + N + 1^ X 

/ /3+7 N+l/ \ \ 

X HcunCi a,l3-i-N-{-l,-y+N + l, -a—^,riA — ( 7V+l-a-|-/3-|-7j , z I 

for all a, iBj'yjrj £ C, (when /3 is not a negative integer) and for any fixed N G X, N > 0. 
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4 A new derivation the confluent Heun's polynomials 



Now we are prepared for an alternative derivation of the confluent Heun's polynomials without using the three- 
terms recurrence relation (|1.3p . i.e., directly from confluent Heun's equation. Indeed, posing the requirement 

Tw(a,/?,7,77) =0, (4.1) 

which is a new form of condition (jl.5bp . we obtain ^^w+i HeunCjv(Q;, (3, 7, 77, z) = 0. Thus under condition (|4.ip 
the (5 jv-confluent- Heun's- function HeunCAf (a, /3, 7, ry, z) becomes a polynomial of degree 

The explicit form (|A.6p of the polynomial condition (|4.ip is given in the Appendix. As seen, it presents an 
algebraic equation of degree [N +1) for the spectral parameter /x of the confluent Heun's equation. Then from 
equation (jA.6P we obtain {N -f l)-in-number roots ^fe=i,...,Af-i-i(a, /9, 7), which yield {N l)-in-number values 
??fe=i,...,Ar+i(a, /3, 7) = ^(q; — /3 — 7 + a/3 — 13^) — /Xjt=i,...,Ar+i(a, 7) of the parameter 77. Hence, the condition 
(|4.ip defines (A^ + l)-in-number polynomial solutions 

PHeunCAr,fc(Q!,/3,7,z) = HeunCAr(a,/3, 7,77fe,2;), /c = 1, . . . , A^ + 1; (4.2) 

to the confluent Heun's equation, each being polynomial of degree N of the variable z. 

Further information about the mathematical properties of the confluent Heun's polynomials 
PHcunC Ar_fc(a, /3, 7, z) one can find in [5]. As a by-product we obtain that under conditions (jl.Sp the con- 
fluent Heun's operator Da,i3.-y,5,ri (|2.ip becomes a quasi- solvable one in the sense of pp] . 

The associate ^Ar-confluent function HeunC^(Q;, 7, 77, z) does not become a polynomial under condition 
(|4.ip . Indeed, the corresponding (^Ar-condition is never fulfilled for it when A^>0 - see the end of the previous 
section. This proofs that the function HeunC* (a, /?, 7, 77, z) is not a polynomial when the ^Ar-confluent function 
HeunCAr(Q!, (3, 7, 77, z) is. 

5 Some comments and concluding remarks 

Both for the general and for the confluent hypergeometric functions a set of simple and universal representations 
of their repeated derivatives in terms of another hypergeometric functions are well known and widely used, see 
for example [21]. 

At present the corresponding mathematical theory is still not developed enough to have a complete picture 
for this problem in the case of different Heun's functions. For general Heun's functions analogous relations 
seem to exist only for some particular cases, see for example |22| 1201 123j . In the known to us literature one can 
not find similar relations for confluent Heun's functions. In the present article we are filling this gap partially, 
taking into account the specific properties of the confluent Heun's functions. 

The general Heun's equation, written in the universal Fuchsian form: 

H"+(^-^ + ^ + ^^)h'+ "^ff ^ g = 0, 7^+5,-fe,=a,+/3, + l; (5.1) 
\z 2 — 1 z—l/aj z(z — l)(z — 1/a) 

was constructed by Karl Hcun in [1 as a generalization of the standard hypergeometric equation, by adding 
one more regular singular point: z = 1/a G (1, oo). Putting 7^ = /3+1, 6^ =7+1, = —a/a, a^P^ = —{n+v) j a, 
q — — /i/a, and taking the limit a — > 0, we obtain the confluent Heun's equation ()1.2p by coalescence of the 
regular singular points z = 1/a and z = 00 in the equation ()5.ip . 

Unfortunately, such coalescence process in the very solutions and corresponding relations is much more 
complicated and, as a rule, not possible. We shall illustrate this fact using the approach of the articles pUlf^[21] . 
Translating the four different regular singular points z = 0, 1, 1/a, 00 of the equation (jS.ip in the compactified 

^Note that the condition 1)4. Ijl coincides with condition l)1.5b| l and its explicit form l)A.6|l up to a nonzero numerical factor. 
It can be shown that the constant 1 13. 8I I is simply related with Starobinsky constant (See 1151 and I17l -|19l.). The equation 114.11 1 
recovers the mathematical meaning of the zero-Starobinsky-constant condition. 
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complex plane Cz to the places of the vortexes a; = 0, 1, 1 + r, t (r being pure imaginary) of the fundamental 
rectangle of periods in the complex plane C^; by the transformation 



p{x) being Weierstrass elliptic function with half-periods ujq = Q,uji = 1/2, (x)2 = 1/2 + t/2,uj3 = t/2, and using 
in addition the substitution 

V-Ca;) =z('"+i)/2(^_l)('i+i)/2(i_a^)(i2+i)/2^(^)^ Zo=7G-3/2,Zi = (5^-3/2,;2 = eG~3/2; (5.3) 
one can transform the general Hcun's equation (|5.ip into elliptic form, identical to the BCl Inozemtsev model: 

1=0 

Here I3 and E are properly chosen constants. This equation is instrumental for deriving the results in the 
articles [501 HSl [H] , where slightly different notations are in use (See, too, the references therein.). Obviously, 
the equation Ij5.4p has no relevant limit a — > 0, since lima^o 1^21 = 00. 

A more fundamental geometrical obstacle to relate the results for general Heun's equation and the results for 
confluent one is that according to the formulae e^ — —j^ei, e2~ ^i+a ^i o^^^ obtains lima^o 62 = 61. Hence, an 
elliptic representation of the confluent Heun's equation does not exists: As a result of the coalescence process, in 
it we have only three different singular points. An analytic map like (|5.2p of the triangle of the singular points 
of the confluent Heun's equation onto a rectangle of the periods of elliptic functions is impossible. Hence, there 
exist essential differences between the properties of general and confluent Heun's functions. Note that some of 
the properties can be reformulated properly after the coalescence. For example, the factor (1 — 02)*^'^+^)/^ has 
a confluent limit exp(az/2). As a result a correspondingly modified substitution (|5.3I) exists. It transforms the 
confluent Heun's equation to Schrodinger-like (non-elliptic) form [3]-[7]. 

A direct consequence of the above consideration is the need to derive most of the basic properties of confluent 
Heun's functions independently of the properties of general Heun's functions. Sometimes one can use similar 
general methods in both cases. An example is the generalized Darboux transformation [20j . used in Section 3 
in a specific way, and invented originally for the equation (j5.4p . 

One can hope that the new mathematical results for the confluent Heun's functions, derived in the present 
article, will be instrumental in different analytical and numerical applications, and especially in the relativistic 
theory of gravity. 
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Appendix 

The confiuent Heun's equation (12. 2[) can be rewritten in the transparently self-adjoint form: 

e-'-'z-P{z ~ lyi- ( e'^'z'+^iz ~ iy+'<^B^) +a(- + ^ + 1^ zi/(z) = /.i/(z). (A.5) 
dz \ dz J \a 2 J 
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Besides, it shows that the natural spectral parameter is /i. Correspondingly, we represent the left hand side 
A7v+i(/i) of condition (jl.Sbp in form of the specific three-diagonal determinant: 



Na fi-q2 + la 2(2+/3) 
{N-l)a fi-q3 + 2a 







fi-qN-i + {N~-2)a 
2a 








{N~l){N-l+/3) 
^i-qN + {N-l)a 
la 









N{N+(3) 
fi-qN+i+Na 



(A.6) 



It turns to be useful in the real calculations [7j. Here q„ = (n — l)(n + /3 + 7). 

A similar representation of the second polynomial condition (jl.Sbp in determinant form was derived in [5]. 

Note that we can develop an alternative consideration, interchanging the places of the regular singular 
points z = and z = 1. Then because of the obvious symmetry of the equation under the change 

{a, 13, 7, S, ^, u, z} — > {a, 7, /3, S, i>,^,z — 1}, the spectral parameter will be —v. 
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